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Abstract 

Theoretical description of off-specular grazing-incidence nuclear resonant scattering of syn- 
chrotron radiation (Synchrotron Mossbauer Reflectometry, SMR) is presented. The recently de- 
veloped SMR, similarly to polarized neutron reflectometry (PNR), is an analytical tool for the 
determination of isotopic and magnetic structure of thin films and multilayers. It combines the 
sensitivity of Mossbauer spectroscopy to hyperfine interactions and the depth selectivity of x-ray 
reflectometry. Specular reflection provides information on the depth profile, while off-specular 
scattering on the lateral structure of scattering layers. Off-specular SMR and PNR intensity for- 
mulae of a rather general multilayer with difi'crent domains, based on a Distorted Incident- Wave 
Approximation (DIWA) are presented. The Distorted- Wave Born Approximation (DWBA) results 
are given in an Appendix. Physical and numerical implications, why using DIWA, are explained. 
The temporal character of SMR imposes specific differences between SMR and PNR. In order to 
reveal the limits of DIWA and to compare the two analytical methods, two-dimensional diffuse 
SMR and PNR maps of an antifcrromagnetic multilayer are calculated and critically compared. 
Experimental 'w — 20' SMR map of a periodic [Fe/Cr]2o multilayer is presented and compared with 
simulations by the present theory. 

PACS numbers: PACS: 42.25-p, 61.10.Kw, 61.12.Ha, 75.25.-|-z 

Keywords: Reflectometry, grazing incidence, diffuse scattering, distorted wave approximation, nuclear reso- 
nance scattering, neutron scattering, multilayers 
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I. INTRODUCTION 



Grazing-incidence reflection from fiat surfaces of x-rays/"^ neutrons^"^ and of soft x-rays 
near tlie absorption edges'''®'^ (x-ray, neutron and resonant magnetic x-ray reflectometry, 
respectively) have been widely used to investigate the chemical, isotopic and magnetic struc- 
ture of thin films and multilayers (ML).^'^"^^ Making use of the sensitivity of nuclear resonant 
scattering (NRS) of synchrotron radiation (SR) to hyperfine interactions, another analyti- 
cal method has been developed that combines the sensitivity of Mossbauer spectroscopy to 
hyperfine interactions with the depth information yielded by reflectometry. NRS of SR is 
performed in the time domain. All nuclear resonances (of typically a few neV linewidth) 
are excited simultaneously by a SR pulse of a few meV bandwidth. Time-differential NRS 
of SR contains the hyperfine interaction information in the quantum-beat pattern of the 
time response that follows the excitation of the system by the synchrotron pulse. Counting, 
however, all delayed (i.e., nuclear resonant scattered) photons as a function of the grazing 
angle of incidence is a special k;ind of x-ray reflectometry that wc shall call henceforth Syn- 
chrotron Mossbauer Reflectometry (SMR).^^'^^ This method is similar to polarized neutron 
reflectometry (PNR) and yields integral hyperfine depth profile and superstructure informa- 
tion. Specular SMR has by now become an estabhshed technique. ^^"^^ 

The specularly reflected radiation from a stratified system is insensitive to the lateral 
structure; it depends solely on the lateral averages of the material parameters^ for the 
coherence volume. If one uses radiation of infinite coherence length, the contrast of laterally 
compensated domains (e.g. the magnetic contrast of an antiferromagnetic ML stack) would 
be absent from the specular reflectivity. However, this is not the case, since the coherent 
averaging has to be performed for the finite coherence volume^^ of the radiation determined 
by the experimental setup and such contributions are to be added incoherently.^^ 

For studying lateral inhomogeneities, such as structural roughness, magnetic domains, 
etc., diffuse scattering, i.e. off-specular rcffcctomctry is used. ('Diffuse' and 'off-specular' 
will be used as synonyms throughout this present paper.) Off-specular non-polarized^^ and 
polarized neutron refiectometry,^'^^ soft-x-ray resonant magnetic diffuse scattering^° and, 
recently, off-specular SMR^° have been used to estimate the domain-size distribution and to 
follow domain transformations in antiferromagnetically (AF)-coupled magnetic MLs. 

The optical theory of specular SMR has been pubhshed using several, somewhat different 
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approaches, ^"'^^"^^ however, to our knowledge, the present attempt is the first one to describe 
nuclear resonant diffuse scattering. Distorted- Wave Born Approximation^^ (DWBA) was 
used earlier by several authors-^^'^-*^"^^ to describe off-specular scattering of neutrons by a 
rough surface. DWBA perturbatively determines the (off-specular) field^'^^"^^ around the 
specular field, the latter being easily calculated, even for general stratified media, by suitable 
matrix methods. ^^'^^'^^'^^'^^"^^ 

The choice of DWBA to elaborate off-specular SMR data seems, therefore, obvious. 
However, unlike PNR, in the SMR case the scattering potential is strongly energy dependent. 
Moreover, SMR is detected in the time domain. Therefore the existing theory can not be 
directly applied. Since the 2D diffuse maps have to be calculated for typically a thousand 
time- (or energy) channels, therefore, in practice, the anyway long calculation time of the 
diffuse maps by an orthodox DWBA would multiply by a thousand and become intolerable 
long for SMR. 

The aim of the present paper is to present an alternative general evaluation algorithm 
for off-specular SMR. A distorted-wave approximation will be used, which is — except for 
exit angles in the total reflection region — accurate enough and, at the same time, it is 
capable of handling the immense computational problem in a reasonable period of time. 
For reasons explained in Sec. Ill, we shall call this approximation 'Distorted Incident- Wave 
Approximation' (DIWA). 

Starting from Lax' general theory^ and from the common optical formalism of polarized 
neutron and Mossbauer reflectometry,^^ an expression will be deduced for diffuse scattering 
of electromagnetic and/or quantum mechanical particle waves on laterally inhomogeneous 
stratified media. From the point of view of specular refiection, the set of discrete atomic 
scattering centers will be replaced by a homogeneous index of refraction n and the scattering 
problem will be traced back to solving the wave equation (far from Bragg reflections).^ In 
the rest of the paper the grazing-incidence limit is studied, for which the index of refraction 
approximation is justified. ^^'^^'^^"^^ 

The material of the paper is organized as follows. In Sec. II, we present a distorted- 
wave approximation without treating the energy dependence of the problem. The presented 
method is a generalization of Vineyard's approximation'^'^ for the case of MLs having domain- 
like inhomogeneities with absorptive anisotropic index of refractions (scattering potentials) . 
In Sec. Ill time-differential and time-integral diffuse intensities are derived for quite general 
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inter- and intra-layer correlations using a domain correlation matrix. The coherent-field solu- 
tion and the cumulative transmissivity for an arbitrary depth is calculated in Appendices A 
and B, respectively. The DWBA formulae for the delayed refiectivity are given in Appendix 
C. In Sec. IV we provide numerical simulations in the DIWA approach for both PNR and 
SMR. Specular and diffuse PNR and SMR scans and maps are compared and the special 
features of the methods are discussed. In Sec. IV D the developed DIWA SMR theory is 
applied to extract the average antiferromagnetic domain size from experimental SMR uj — 26 
maps of an [Fe/Cr] epitaxial ML. 

II. OFF-SPECULAR SCATTERING 

The model systems of the present paper are MLs and thin films, having lateral inhomo- 
geneities on the mesoscopic scale. Unlike the case of surface roughness, these inhomogeneities 
(domains) will be assumed to be much larger than the atomic distances. In each homoge- 
neous part around position r an index of refraction n (r) is defined. Since the elements 
of n for both slow neutrons and X-rays differ only shghtly (typically 10~^ — 10~^) from 
that of the 2x2 unit matrix /, the small parameter, namely the generalized susceptibility, 
X (r) = 2 [n (r) — /] can be conveniently defined. 

Using the index-of-refraction approximation, in each homogeneous part of the system, 
the solution of the homogeneous wave equation 

[A + A;'/] * (r) = -ifc^x (r) * (r) (2.1) 

yields \[' (r), representing the two components of the photon field or the neutron quantum 
mechanical spinor state at position r, with k being the wave number in vacuum, x is 
simply related to the coherent forward-scattering amplitude / by x = ^p^/ where N is the 
number of scattering centers per unit volume.^^'^^ For photons / = /e + /n is the sum of the 
electronic and nuclear scattering amplitudes.^^ For neutrons / = /nuc + /magn is the sum of 
the nuclear and magnetic scattering lengths. In Eq. (2.1) k'^x plays the role of the optical 
scattering potential (scattering- length density in the neutron literature). For photons x is 
the susceptibility and, for sake of simplicity, throughout the paper we shall use this term in 
its general sense. 

Describing the system as a stack of (possibly laterally inhomogeneous) layers, we compose 
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the susceptibility by 

X{r) = f2^i{r\\), (2.2) 
1=1 

as the sum of the susceptibihty functions of the individual layers I {I — 1..S, the last layer 
S being the substrate) depending solely on the in- plane coordinate ry . 

If the homogeneous parts of the system are large compared to the wavelength wc may 
assume that the exact solution "if (r) is close to the solution ^'coh (r) of the coherent (specular) 
field equation^ 

s 

[A + k'l] *,oh (r) = -k' J2 ^^*coh (r) , (2.3) 

1=1 

which is obtained from Eq. (2.1) by replacing the susceptibihties xi by the average 
susceptibility Xi of each layer I. In order to deduce a perturbative equation, the sum 
—k'^ ^X;\E' (r) is added and subtracted on the right-hand side of Eq. (2.1) 

s s 

[A + k'l] ^ (r) = -k'J2xi^ (r) -k'Yl M ' ^i] ^ (r) • (2-4) 

/=i 1=1 

For homogeneous layers Xi (ry) = Xi ! ^•^•) the second sum vanishes on the right, so that Eq. 
(2.4) reduces to Eq. (2.3), the basic equation of specular refiectometry.^^'^'''^^''^^ 
The general solutions of Eq. (2.4) are searched for in a form 

*(r) = *coh(r) + *off(r) (2.5) 

where ^coh (r) is the coherent field, which vanishes in any non-specular direction, and ^ofr (r) 
is the off-specular field. Substituting Eq. (2.5) into Eq. (2.4) and taking into account Eq. 
(2.3) we obtain 

[A + k'l] (r) = -k' [xi (rii) - xi] *coh (r) - k'J^Xi (r||) *off (r) . (2.6) 

1=1 1=1 

The coherent field ^^coh (r), solution of Eq. (2.3), is obtained by the optical method^^ as 

*eoh (k, r) = T {k±, rx) *^^exp (zk||r||) (2.7) 

(see Appendix A) where _L and || denote the plane-perpendicular and in- plane components 
of the respective vectors and ^™ is the amplitude of the incident plane wave of wave vector 
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k = {k±, k||). Here we introduced the 'cumulative transmittance' of the reflecting fllm from 
surface to a depth r± by: 

T {k^, r^) = LPi] {k^, r^) [I - R^^ (k^)] + L^^l (A;^, r^) [/ + R^^ (k^)] . (2.8) 

Rsp{k±) is the 2x2 specular reflectivity matrix of the system^-*^, L^'^^^ {k±,r±) and 
L'^^' {k±,r±) are the respective 2x2 submatrices of the 4x4 characteristic matrix,^-^'^^ 
L at depth r± for an incoming plane wave defined by k. 

The physical interpretation of the inhomogeneous wave equation (2.6) is seen from its 
right-hand side where the first term gives the source of the off-specular radiation and the 
second term describes the off-specular field scattered by the entire ML. The off-specular field 
arises from the coherent field at the lateral inhomogeneities, i. e. from the regions where the 
susceptibility x differs from its average value x- Eq. (2.6) can be solved iteratively to the 
prescribed accuracy. 

As a first approximation, the second term of the right-hand side of Eq. (2.6) is neglected 
so that the solution can be obtained using the Green-function technique, 

*off (k, r) = ^ E / ^^^^1^ [Xi (rli) - Xi] ^coh (k, r') , (2.9) 

where R— |r — r'| . The approximation requires ||^coh (r)!! ^ H^ofr which condition is 
fulfilled for magnetic MLs of large enough homogeneous domain size in the vicinity of the 
specular direction^*^'^^ so that the exact solution ^ (r) is close to the coherent field \E'coh (r)- 
When neglecting the second term in Eq. (2.6) the scattering of the off-specular field is 
neglected. Therefore the present distorted- wave approximation breaks down for exit angles 
near the critical angle of total refiection. 

Far from the scatterer, the Praunhofer approximation 

'-^^^^'-^^e.p{-z].'r') (2.10) 

is applied with k' being the wave number vector of the emerging plane wave, with which the 
final expression for the off-spccular field is 

*off (k, r=jr^= exp (ikr) J] Si (K||) 7] {k^, k'^) ¥\ (2.11) 

with K|| being the in-plane component of the momentum transfer vector K = k' — k, and 

Ti{k^,k'^) = -^J dr^cxp{-ik'^r^)T{k^,r^) (2.12) 



is the Fourier integral over the one-dimensional interval Zi of layer I. The expression 

Si (K||) = i- y d^rii exp (-zKy r,,) {xi (ry) - xl] (2.13) 

is the two-dimensional Fourier transform of Xi {^\) ~ X- Ti {k±, k'j_) that can be analytically 
calculated (See Appendix B). 

The off-specular scattered intensity /off = (*off, ^off) is 

^off = E (*'^ {k^, k'^) Sj (K||) S, (K||) T, {k^, k'^) ^>'-) , (2.14) 

w 

which, for arbitrary incident polarization, can be rewritten as 

where p is the polarization density matrix of the incident radiation.^^ From the convolution 
theorem, it follows that the Fourier transform Cw (Ry) of 

C„.(Ky) = (27r)5'/(Ky)5;.(Ky) (2.16) 

is the cross-correlation function of the susceptibilities between layers / and 

(Ry) = J d^ry [xi (R|| +ry) -Xi]^ [xf (ry) -Xv] ■ (2-17) 
The final result then becomes 

a convenient expression for randomly distributed lateral inhomogeneities. We note that 
the off-specular intensity Iqs — /off (Ky, k±, k'^ is a function of Ky, k±_ and A;^, a notation 
dropped in the calculations. The corresponding values of Ky, k^ and k'^ can be given for the 
chosen experimental geometry. The diffuse intensity expression in Eq. (2.18) differs from 
that of DWBA,^° an issue to be discussed in Sec. Ill and in Appendix C. 

A possible experimental realization of the off-specular reflectometry is the so-called 'o;- 
scan' geometry where the detector position is set to 29 and the sample orientation uj on the 
goniometer is varied with the sample normal remaining in the scattering plane. For a;-scans, 
the in-plane components of the momentum transfer vector Kn and the plane-perpendicular 
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component of the wave vector of the emerging wave /c^can be expressed by the grazing angle 
9 and by uj: 

K\\^2ksinesm{uj -9) , (2.19) 

k^^ksinuj, (2.20) 

fc^ = -fcsin(2^-a;). (2.21) 

One may observe that K\\ vanishes at the specular condition uj — 6. A two-dimensional 
representation of the full diffuse scatter is the 'a; — 2^'-map. A further, widely used, ar- 
rangement is the so-called 'detector scan' geometry in which the angle of incidence is 
fixed and the scattered intensity is recorded as a function of the exit angle ^out- In case of 
a detector scan, Eqs. (2.19) and (2.24) read 

A'll = k (cos 6'out — cos 6'in) , (2.22) 
A;sin^in, (2.23) 
-fcsin^out, (2.24) 
An alternative representation of the full diffuse 



k^^ 

with the specular condition being ^out = ^in- 
scatter is the '^^m— 6'out' niap. 



III. OFF-SPECULAR SMR 



A. Time-differential scattered intensity 

Eq. (2.18) for the off-specular intensity is equally valid for neutron and x-ray scattering 
where the time of the scattering process is negligible. However, in case of nuclear resonant 
scattering we detect the time response 

oo 

*off {r,t) = -4= / dE *off {r,E) exp {-iEt/h) (3.1) 

— oo 

after the synchrotron pulse, "'"^ which is the Fourier transform of the energy-dependent ofT- 
specular field. Close to a Mossbauer resonance, both the susceptibilities xi (^\\ > ~ Xi i^) 
and the coherent field ^^coh (i"', E) are strongly energy-dependent.^^ Therefore, through Eqs. 
(2.8), (2.12) and (2.13), the Si (K||,£;) and the Ti{k±,k'^,E) quantities carry an energy 
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dependence, too. Consequently, Eq. (2.18) can no longer be applied to calculate the off- 
specular intensity. 

A possible workaround of this problem is to define a distribution function flf (ry) of 
homogeneous regions of type /i = 1, .., M of layer /. This function (ry) characterizes the 
homogeneous regions of layer / of an energy-dependent susceptibility (E). Over the region 
of type fj, in layer I the distribution function flf (ry) — 1 otherwise Of (ry) — 0. Any point 
ry along the surface of layer I is related to one domain type, therefore 

E^'h) = l- (3-2) 
The total inhomogeneous susceptibility is the sum 

M S 

x(r,^) = EE^r(r||)x''(^) (3.3) 

where the space- and energy-dependent terms in x {'"^^E) have been separated. The average 
susceptibility within layer I is 

M 

Xi{E) = J2vh'iE) (3.4) 
where the fractional domain area is 

77f = Af/A (3.5) 

i.e., the ratio of = / cPr\i Of (ry) , the total area of the homogeneous part of type within 
layer Z, and A, the area of the ML. Since the domains fully cover the layers, the condition 

E^r = l (3-6) 

is fulfilled. Using (3.3) and (3.4), Eq. (2.13) becomes 

M 

with 

Wl" (Ky) = i- 1 d^rii exp (-^Kyr,,) [Of (r,,) - r^f] (3.8) 
and, finally, the energy-dependent off-specular field in the Fraunhofer approximation is 

^off {r,E) = J^y exp (tkr) E (K||) (E) (E) (3.9) 
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while the off-specular intensity is 




(K||)Tr {E)' {E) p 



(3.10) 



The matrices for the homogeneous region /x, / 



rf {E)^x''{E)Ti{E) 



(3.11) 



are the products of the homogeneous solution T; [E) and the susceptibility {E) of that 
region. For the sake of brevity, henceforth the dependence of 7] and on both k±_ and 
/c^will not be explicitely written. 

As we already indicated in the previous section, Eq. (2.18) and, consequently, Eq. (3.10) 
differ from the usual expression of diffuse intensity of DWBA as published earher in the 
literaturc^^'^^"^^ in such a way, that the "distortions" are only considered on the incident 
path before scattering, the exit path is left undistorted. The DWBA expression is given and 
discussed in Appendix C. Physically speaking, the difference between the two approxima- 
tions is that DIWA only takes distortions into account on the incident wave, while DWBA 
on both incident and emerging waves. Therefore a DIWA calculation is considerable faster 
but it is not expected to be invariant with respect to exchanging the source and detector 
(a condition widely called 'reciprocity'"'^"'^^). Should, nevertheless, reciprocity be physically 
justified under some conditions, we may take advantage of calculating only half of the '^^m— 
^out' and the 'a; — 2^' maps along with mirroring one side of the = ^out line and the 
uj — 9 lines onto the other, respectively. We shall see that the DIWA calculation is quite 
accurate at one side of these lines, a fact finally resulting both in saving another 50 % of 
the computation time and in a high accuracy of the calculation of the full diffuse scatter. 

The geometrical (or 'domain') correlation function C^'^f between layers / and /'and ho- 
mogeneous parts n and p'oi the layers with definition 



may describe quite general structural and magnetic intra- and interlayer correlations like 
correlated layer growth, or magnetic-magnetic correlation (e.g. antiferromagnetic multilayer 
domains) and even magnetic-structural correlations. Similarly to Eq. (2.17), the direct 
space correlation function can be written as 



C/;f'(K||)=(27r)Ty,-(K||)* < (K||) 



(3.12) 




(3.13) 
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Notice that both the correlation function and its Fourier transform are symmetric with 
respect to the simultaneous exchange of the domain and layer indices: 

C (R||) = (R||) • (3-14) 
An important consequence of Eqs. (3.2) and (3.6) is 

(R||) = ' (I^ll) = 0- (3-15) 
Using Eq. (3.13) the correlation function at R|| = reads 

C' (0) = r^!;^' - «' (3.16) 

where 



\ [ d\||Qf (r||)Qj;'(r||) (3.17) 



Vw ^ 

is the fractional overlap of homogeneous parts of types n and fi' in layers / and /', respectively. 

In Eq. (3.10) the geometrical correlation is separated from the energy-dependence and, 
therefore, it can be applied for time and energy domain experiments alike. The Fourier 
transformation can be performed so that, using Eqs. (3.1) and (3.10), the time-dependent 
intensity becomes 

^off (^) = 1^ E Tr k {ty G^: it) p] (3.18) 



4^2 



where 



oo 

it) = ^ J dETf {E) exp i-iEt/h) . (3.19) 



-oo 



B. Time-integral scattered intensity 

The time-integrated intensity is /"^' = j dt los (t) where ti and t2 define the starting and 

ti 

finishing time of the time window of the counting after the synchrotron pulse. Applying 
Eqs. (3.18) and (3.19) 

oo 

= E / E (K||) Tr [rf {E + may {E) p] (3.20) 

m=-oo _^ ll'iifj,' 

where e = ^ with tbunch being the time interval between the synchrotron bunches, h 

I'bunch 

is the Planck constant, Sm is the m}^ discrete Fourier component of the periodical time 

oo ^ X 

window function S (t) = ^ s^exp [i^^tj oi the experiment defined by S (t) — 1 for 

m=— oo ^ ' 

^1 <t <t2, otherwise S {t) after each synchrotron bunch. 
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IV. MODEL CALCULATIONS AND COMPARISON WITH EXPERIMENT 



Applying the above theory, off-specular time-integrated SMR curves and maps were 
calculated and, on one hand, compared with PNR curves and maps simulated with the 
same DIWA theory and, on the other hand, compared with experimental off-specular 
SMR data. In order to treat realistic problems, at this point further specification 
of the studied system and the experimental conditions will be undertaken — without 
restricting the generality of the discussion. The actual calculations were performed for the 
MgO/ pPe (2.62 nm) /Cr (1.28 nm)]2Q multilayer structure. No interface layers and no rough- 
ness contribution were considered. Although, as mentioned, the Cjti^ functions may describe 
a variety of domain correlations, simulations were performed for the case of strongly-coupled 
layer antiferromagnet in remanence. By layer antiferromagnet we mean even number of 
magnetic layers of identical thickness resulting in zero net magnetization of the ML stack 
in remanence. Except for the trivial case of full in-plane saturation, the magnetic layers in 
an AF multilayer stack are broken into domains of different orientations. In remanence, 
we assume 180° domain walls of negligible thickness, i.e. the sublayer magnetizations also 
vanish (for further details see below). Strong coupling implies a strict plane-perpendicular 
domain correlation throughout the ML stack, a structure, in which the top layer unequivo- 
cally identifies the domain structure in the lower ones (say '+' and ' type domains). The 
actual functional form of such AF domain correlation is given in the next section. 

For the SMR simulations, the scattering geometry was selected so that the layer mag- 
netizations lay both in the plane of the film and of scattering, i.e. parallel/ antiparallel 
to the in- plane component of the wave vector, ky, a condition for the appearance of the 
SMR specular AF reflection. All SMR curves and maps were calculated for the 14.4-keV 
Mossbaucr resonance of ^''Fc (A = 0.086 nm) for hyperfine magnetic fields of ±33.08 T for 
the'-|-' and '— ' type domains, respectively. The isomer shift and quadrupole sphtting were 
set to zero (parameters for bcc iron at room temperature). The electronic susceptibilities 
for the A = 0.086 nm x-rays were taken from the Berkeley web site^^ for the various ele- 
ments. The synchrotron bunch time ^bunch, and integration boundaries ti and t2 were chosen 
according to the actual experimental values (see below). 

For PNR simulations, monochromatic neutrons of wavelength A = 0.4 nm were assumed 
and the layer magnetizations were set parallel/antiparallel to the neutron spin and per- 
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pendicular to k||. The nuclear scattering lengths used for Cr, ^"Fe and MgO are 3.6, 
2.3, and 11.2 fm, respectively.^^ The magnetic scattering length bjn{z) — Cijl{z), where 
C = ro7/2 = 2.69542 fm//XB, where Tq is the classical electron radius, 7 = 1.91304 is the 
magnetic moment of neutron in nuclear magnetons, ^{z) is the average magnetic moment 
per atom per unit volume at depth z, ji-Q is the Bohr magneton. Numerically, hy^{z) was 
set to ±5.93 fm for the '+' and '— ' type domains within the ^^Fe layers and zero otherwise. 

The theory presented above was implemented in, and specular and off-specular intensity 
curves and maps (both SMR and PNR) were simulated by the data evaluation computer 
program EFFI (Environment For Fitting), which is freely downloadable.^^'^° 

A. Consequences of the finite coherence length 

In the case of a 1:1 surface coverage of the '+' and '— ' type domains, laterally averaging 
the magnetizations within a layer — as mentioned in the introduction — leads to a loss 
of the AF contrast and no AF Bragg peak in the specular reflectivity appears. However, 
when the lateral size of the domains is bigger than or comparable to the lateral coherence 
length^^ of the applied radiation, a net layer magnetization is sampled within the coherence 
area and such intensity contributions are to be added incoherently. ^^'^^'^^ Consequently, a 
magnetic contrast appears in the specular reflectivity. In order to account for the effects of 
the finite coherence volume, the domain size is to be related to the lateral coherence area, 
the projection of the coherence volume^^ to the top magnetic layer. (The lateral coherence 
length was reported to be 0.1 to 30 /im for neutrons^^ and similar values can be derived for 
nuclear resonant photons. 

In order to account for the effects of the finite coherence volume, the domain size is to 
be related to the coherence area. In order to account for the effect of the finite coherence 
length we redefine the fractional area of a domain type, r]i that was defined above for incident 
radiation with infinite coherence lengths. Indeed, as a first approximation, one can use Eq. 
(3.5) inside the coherence area within the top layer by exchanging A, the total area of the 
film, for the coherence area and relating the surface of type domains to the coherence 
area. In case of an antiferro magnet the domain type index yU identifies the '+' and ' — ' type 
domains (/x = +, — ). Using Eq. (3.6) one can introduce a specific magnetic bias parameter 
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77 for the layers of even and odd index by the definition 



V = Vli+i = 1 - V2i+i = % = 1 - Vli 



(4.1) 



with i being an integer. 

In Fig. 1 simulated specular SMR intensity curves are displayed for different values of rj. 
As expected, the r^-dependence is restricted to the (= ^^out) regions of 1/2- and 3/2-order 
AF Bragg reflections. As explained, there are no AF peaks for rj = 0.5 and the magnetic 
contrast increases with increasing dominance of either '+' or '— ' type domains. Since the 
coherence area is smaller then the illuminated area and it may have an arbitrary position, 
which can be identified by r), the integration in Eq. (2.17) is performed as a function of rj. 
The off-specular intensity in (2.18) is calculated as 



where p {rj) is the normalized probability density of having a (magnetic) domain bias of rj. 
We note that the probability function p [r]) depends both on the lateral domain structure and 
on the measurement geometry (source and slit dimensions, etc.). For the sake of simplicity, 
we consider here the case when the probability density p [r]) sharply peaks near and 1 
since strong specular SMR AF Bragg peaks can indeed be observed experimentally. Since 
the NRS spectra of a magnetic layer with /c-parallel or a /c-antiparallel directions of the 
hyperfine field can not be distinguished, the SMR AF Bragg peak is of the same shape and 
intensity irrespective of whether the AF structure starts with a A;-parallel or a A;-antiparallel 
layer on the top of the stack. Therefore the curves in Fig. 1 are identical for domain bias 
Tj and 1 — rj. 

The specular curves in Fig. 1 are typical for an AF multilayer. For the above model struc- 
ture and wavelength corresponding to the ^^Fe Mossbauer resonance energy the structural 
and AF Bragg peaks show up at = 11.3 and = 6.7mrad, respectively. The Kiessig 
fringes, characteristic of the total thickness of the ML, appear in all specular curves, how- 
ever, due to the strong nuclear resonant absorption in the "''^Fc layers, their amplitude is 
strongly attenuated. Below the critical angle, the delayed reflected SMR intensity tends to 
zero at zero incident angle, ^^'^^ whereas a non-resonant x-ray or PNR curve behaves 'nor- 
mally' (i.e. tends to unity, corresponding to total reflection). This difference, as we shall 
see, leads to an augmented effect on the diffuse intensities. 



1 




(4.2) 
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B. Lateral correlation: domains 



The domain structure within the top layer of the AF multilayer stack will be assumed 
laterally isotropic and statistically characterized by the correlation function of the top layer 
magnetization. For simplicity, we assume an exponential function and interpret the average 
domain size as the correlation length ^ of the exponential correlation function^° 



'^tl^' (ril) = (-1)'+''+^-'-' 77 (1 - ,7) exp [-'j^ 



(4.3) 



where 5 is the Kronecker delta symbol. We assume the same correlation function for each val- 
ues of indices observing an alternating sign corresponding to the strict plane-perpendicular 
AF correlation and an in-plane '+'/'—' domain model. The Fourier transform of Eq. (4.3) 
reads 

(K||) = (-1)'+^'+-^..' r^il-r^) (4.4) 

which is then substituted into Eqs. (3.10) and (3.20) for calculating the off-specular PNR 
and SMR intensities, respectively. 



C. Off-specular scattering: PNR and SMR 



Figs. 2 and 3 show simulated a;-scans for the above model multilayer at 29 — 52 mrad 
and 29 = 13.4 mrad corresponding to the 1/2-order AF Bragg angle for different domain 
correlation lengths ^ using a single domain bias parameter 7] = 0.1, for PNR and SMR, 
respectively. Applying Eqs. (C6)-(C9) the kinematical approximation (first Born approxi- 
mation, BA) is also shown. The curves were symmetrized by substituting the mirror image 
of a; < ^ in place of the u > 9 region. The grounds of such substitution will be explaned be- 
low. The BA curves display the same shape in the entire ^-range and the off-specular scatter 
width is inversely proportional to ^. The BA and symmetrized DIWA (called, henceforth, 
'Symmetrized Distorted Incident-Wave Approximation", SDIWA) cm^vcs for PNR overlap 
almost in the entire angular range, the differences only develop near the critical incident 
and exit angles, i.e. at the Yoneda wings. For SMR, the SDIWA curves change their shape 
for different correlation lengths, since, as expected, the plane-perpendicular electronic and 
hyperfine depth profile, the energy- and polarization-dependent absorption embedded in the 
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rf functions have a strong influence on the shape of the off-specular scatter. Unlike in case 
of PNR, therefore, identifying the correlation length with the inverse width of the diffuse 
SMR scatter - without a proper evaluation of the entire reflectivity curve - can not be justi- 
fied. Note, that the simulations were performed for the experimentally feasible neutron and 
nuclear photon wavelengths of A„ = 0.4 nm and Asrpe = 0.086 nm, respectively. Although 
the scattering amplitudes are comparable, the angular range between the Yoneda wings, 
consequently the angular width of validity of the BA approximation is considerably wider 
for the longer wavelength neutrons. This, however, does not effect the above statements on 
the quahtative differences between the PNR and SMR maps. 

Figs. 4b and 5b display simulated two-dimensional '^^in — ^^out' niaps for ^ = 1 /xm corre- 
lation length with a single domain bias parameter = 0.1 in the same DW approximation, 
corresponding to Eqs. 3.10 and 3.20 for PNR and SMR, respectively. The intensity is max- 
imum along the diagonal specular line, and a broad diffuse intensity is observed around the 
half-order Bragg peaks. Similarly to the a;-scans, the < ^out semi-plane was mirrored 
onto the > ^out semi-plane. Only the diffuse SDIWA and BA maps are displayed for 
PNR in Fig. 4. /"'""and /""*" maps are not shown since the spin-flip scattering vanishes due 
to the '+'/'—' domains being assumed parallel/antiparallel with the incident neutron spin. 
Moreover, due to the small absorption of the neutrons, the -|— I- and maps (correspond- 
ing to the same layer structure with reversed layer sequence), are practically identical (and 
therefore not shown) . It is not surprising that the Kiessig fringes are observed in both PNR 
and SMR diffuse scatter since the source for the diffuse intensity is the specular field. Due 
to the negligible absorption of the neutrons, the Kiessig contrast is stronger in PNR than 
in SMR. We may observe further symmetries in the PNR curve. As expected, the diffuse 
intensity around the structural Bragg node at 9in = 9out — 52 mrad is missing, since the 
diffuse scattering here is purely magnetic origin and the magnetic contributions cancel each 
other at the momentum transfer corresponding to the first order Bragg peak. Moreover, 
for ^'^Fe/Cr for a 2:1 layer thickness ratio the 3/2-order Bragg refiection is forbidden, conse- 
quently, the node for the 2.62 nm(^^Fe)/1.28 nm(Cr) layer thickness ratio in Figs. 4, and 5 
are very weak. 

The ^''Fc ^^in — ^^out SMR map in Fig. 5b shows a number of features different from 
that of PNR. Due to the dominance of energy-dependent absorption for resonant x-rays the 
SMR BA and SDIWA maps drastically differ. E.g., weak structural Bragg wings around 
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^in = ^out = 11.25 mrad appears, but those fall below the typical experimental background 
level. Below the critical angle (unlike in cases of PNR or non-resonant x-rays) , the specular 
SMR intensity, the source of the diffuse scatter, rapidly decreases to zero (c/. Fig. 1). Since, 
independently, due to the stronger absorption, the Kiessig fringes arc suppressed in the SMR 
map in Fig. (5b), the intensity does not oscillate near Yoneda wings, which, as a result, 
appear broadened as compared to the PNR pattern in Fig. (4b). As a consequence, the 
total-reflection peak, being in the critical region, is somewhat difficult to distinguish from 
the 1/2— order AF Bragg node in Fig. (5b). 

So far, by discussion of simulated SMR and PNR-scans and maps, we have shown, that 
the presented DIWA approach satisfactorily describes both the PNR and SMR off-specular 
intensity — except for exit angles below and around the critical angle of total reflection. 
In Appendix C the more exact DWBA SMR intensity formula is also given, which takes 
into account the distortions on both incident and emerging waves. However, the typical 
computation time needed for the DWBA calculation renders its usage completely impracti- 
cal for fitting of diffuse synchrotron Mossbauer reflectograms even on advanced present-day 
computer architectures. Indeed, Eqs. (3.10) and (C6), differ only in the definition of Ff 
according to Eqs. (3.11) and (C7), for DIWA and for DWBA, respectively. Counting the 
number of the 2x2 complex matrix multiplications one can compare the speed of the two 
algorithms. On the one hand, in the DIWA equations (3.11) and (B7), one has to perform 
nine complex 2x2 matrix multiplications. On the other hand, in the DWBA Eq. (C7) 
Np X 13 complex 2x2 matrix multiplication have to be performed. Np is typically ten. 
One may, therefore, conclude that DIWA, as derived here, is, at least, 10 x 13/9 fa 14 
times faster, than DWBA. Since, as previously explained, a single diffuse SMR map requires 
as much more computing time as the number of energy channels in the Mossbauer spec- 
trum (typically 1024), this further factor in computing time would make such calculations 
extremely tedious. Therefore an alternative approach was followed here. 

In case of PNR, as a consequence of the neghgible absorption of neutrons (Hermitian 
scattering potential), reciprocity is fulfilled. Conversely, the resonant Mossbauer medium 
is absorptive and gyrotropic and reciprocity can not be proved to hold in general. ^^'^^ In 
fact, the violation of reciprocity in scattering on absorptive and grotropic media is still the 
subject of both experimental^^ '^^'^'^ and theoretical^^'^^ works. 

Nevertheless, for the special case under discussion, when all hyperfine magnetic fields are 
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aligned parallel/antiparallel to the wave vector, one may verify that reciprocity exactly holds. 
Indeed, in this case, the scattering amplitude is diagonal on the circular polarization basis 
throughout the whole multilayer and, therefore, Eq. (2.1) becomes uncoupled. Consequently 
the conventional proof of reciprocity to both eigenmodes is straightforward"^^'^^ and the '^m — 
^out' SMR maps become symmetrical. In view of the above, we presribe reciprocity for both 
PNR and SMR and mirror the < ^out semi-plane onto the > ^out semi-plane. By doing 
so, the DWBA accuracy of the < ^out semi- plane is achieved by the present (faster) DIWA 
algorithm on the entire '^in— ^out' plane (with a further decrease of computing time by a 
factor of two. This symmetrizing procedure (SDIWA) was used to simulate the DIWA curves 
and maps in this paper (cf. Figs. 2, 3, 4, 5 and 6). 

D. Experimental results and discussion 

Two-dimensional experimental and simulated 'a; — 2^' SMR maps of a [^^Fe/Cr] an- 
tiferromagnetic multilayer are presented in Fig. 6 in the vicinity of the antiferromag- 
netic (1/2— order) Bragg reflection (region marked by dashed lines in Fig. (5)). The 
MgO(001)/p''Fe/Cr]2o ML was prepared by molecular beam epitaxy at the IMBL facility of 
IKS Leuven, Belgium. Preparation and characterization of the sample has been described 
elsewhere. The layering was verified epitaxial and periodic, with thicknesses of 2.6 nm 
and 1.3 nm for the ^^Fe and Cr layers, respectively. SQUID magnetometry showed a satu- 
ration field of 0.9 T and AF coupling between neighboring Fe layers. According to previous 
studies on this multilayer, the Fe magnetizations at remanence align along the (100) 
and (010) perpendicular easy directions corresponding to the respective (110) and (110) 
directions of the MgO substrate. Experimental realization of alignment of domains along 
the fc- vector was achieved in the following way. First a magnetic field of 1.6 T was applied 
perpendicular to k in one of the two equivalent in-plane directions of easy magnetization of 
the MgO(001)/p^Fe/Cr]2o ML, then the field was relaxed to remanence. By this procedure, 
due to the antiferromagnetic coupling between the layers, the sublayer magnetizations be- 
came aligned in the perpendicular easy direction, the two types of AF domains being only 
different in the top-layer magnetization direction. 

Experimental u — 26 SMR maps were recorded at the BL09XU nuclear resonance beam 
line^® of SPring-8, Japan, using the 14.4 keV Mossbauer transition of ^^Fe by performing 
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sequential a;-scans in a 2^-range. The synchrotron was operated in the 203-bunch mode, 
corresponding to a bunch separation time of ^bunch = 23.6 ns. The SR was monochromatized 
by a Si(4 2 2) /Si(12 2 2) double channel-cut high-resolution monochromator with a resolution 
of 6 meV. The specimen was mounted in grazing-incidence geometry with a sample-to- 
dctcctor distance of 46 cm and receiving slit width of 0.1 mm. The delayed radiation was 
detected using three Hammamatsu avalanche photo diodes (APD) behind one another in 
order to increase detecting efficiency. The delayed photons were time integrated using the 
time windows given by ti — 1.97 ns and t2 — 21.63 ns. For these special APDs the dead time 
was as small as below ti. Fig. 6 shows the two-dimensional experimental (a) and simulated 
(b) 'tu — 2^' SMR maps in the vicinity of the 1/2-order Bragg position. In order to avoid 
transformation of the angles varied in the experiment, instead of the 6^^ — 9out plane, data 
and simulations are displayed in the uj — 29 plane. Along the diagonal, the specular line 
was added taking the experimental receiving slit width into account. By trial and error 
and visual comparison, the best correlation function (4.4) was found with the correlation 
length of ^ = 1.0 /im. The simulation is in rather good agreement with the experimental 
data. A detailed discussion of experimental diffuse SMR scans including external field and 
field-history dependence of the domain structure will be reported elsewhere. 



V. SUMMARY 



In summary, expressions for the diffuse scattering intensity of grazing-incidence nuclear 
resonant scattering of synchrotron as well as polarized neutron radiation have been derived 
in a distorted-wave approximation. Distortion only of the incident wave was taken into 
account. In a common optical formalism,^^ grazing-incidence (specular) x-ray, polarized 
neutron, nuclear resonance reflection, and grazing-incidence diffuse scattered intensity were 
calculated in terms of (geometrical) domain correlation functions and the specular field 
depth profile. The formula describes scattering by domains of rather general types with 
intra- and inter-plane correlations and is not restricted to the treatment of random lateral 
roughness. It was shown that, since the off-specular scatter shape is strongly dependent on 
the in-plane correlation length of a single exponential correlation function, without properly 
accounting for the specular fleld depth proflle, no conclusions can be drawn on the shape 
of the domain correlation function. By prescribing the reciprocity (mirror symmetry of 
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the — ^out maps) the inadequacy of DIWA for critical exit angles was ehminated. In 
addition, reciprocity was shown to exacly hold for the most widely used case of SMR, i.e., 
when the layer magnetizatios are parallel/ antiparallel to the photon wave vector. The code 
based on the presented theory is suitable for simulation of diffuse SMR maps in a feasible 
calculation time, which is not yet the case for the DWBA formulae presented in Appendix 
C. Off-specular SMR and PNR a;-scans and — ^out maps of an antiferromagnetic [Fe/Cr] 
multilayer were calculated and compared to each other in order to show the different features 
of diffuse scattering of nuclear resonant synchrotron radiation and polarized neutrons. As 
an application of the presented theory, cu — 29 nuclear resonant diffuse scattering maps of an 
MgO/ pFe (2.62 nm) /Cr (1.28 nm)]2o muhilayer were simulated and, using an exponential 
in-planc correlation function, a layer magnetization correlation length of ^ = 1.0 /^m was 
derived for the ML demagnetized from easy axis saturation to remanence. 

APPENDIX A: GENERAL SOLUTION OF THE COHERENT FIELD EQUA- 
TION 

The solution of Eq. (2.3) was given in Refs.^-'^'^''''^^ where, using the derivative field 
$coh (^±) = i'ik sin 6)^^ "^'coh (^J-)) ^^e second-order differential equation regarding to \&coh (r), 
was replaced by a set of first-order differential equations, providing the solution 



where L is the 4x4 characteristic matrix of the system,^^'^"'^^ k± — ksinO is the plane- 
perpendicular component of the wave number vector of the incident plane wave, which 
latter dependence we drop in this appendix. Here ^ and $ are coherent fields; the notation 
'coh' has been dropped. The physical meaning of Eq. (Al) is that there exists a linear 
connection expressed by the characteristic matrix L between the fields at depth r± = and 
at an arbitrary depth r±. Taking into account the boundary conditions, the field at the top 
surfaces of the system {r± — 0) is 



i.e., the sum of the incident and the reflected -Rgp^™ waves so that Eq. (Al) reads 




(Al) 



^ (0) = + i?,p^i^ 



(A2) 




(A3) 
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where the concept of impedance tensors by Ref. ^ was used, taking into account that the 
fields at r_L = are in vacuum (see Eqs. (21) and (22) of Ref.^^). Expressing the second 
component from Eq. (A3), the field at an arbitrary depth r± we have 

* (r^) = [lP^I (r^) (7 - R,,) + L^^l (r^) (7 + 7t:,p)J (A4) 

and using the notation T (r±) = L^^l {r±) (7 - R^p) + L^^l {r±) (7 + R^p) the solution of the 
three-dimensional homogeneous wave equation is 

*eoh (k, r) = T (k±, r±) exp (ik||r||) . (A5) 

APPENDIX B: CALCULATION OF THE Ti {k^,k'^_) FOURIER INTEGRALS 

In this appendix the analytical calculation of the integral (2.12) is given. The 4x4 
characteristic matrix^'''^^ of an arbitrary homogeneous multilayered film with layers / — 
1, S is 

L = Ls...L2Li (Bl) 

where 

( cosh (kdiFi) lFismh{kdiFi)\ 
Li = \ I (1^2) 

\ xFj^^ sinh {kdiFi) cosh (kdiFi) I 

is the characteristic matrix of the l^^ homogeneous layer^^'^^ with di being the thickness of 
the l^^ layer, x — isinO and the 2x2 matrix Fi — •>/— 7sin^ & — Xi- We note that L depends 
on k± — ksinO, which dependence is not indicated in this Appendix. At depth r± measured 
from the top of the ML the position vector points into layer j < S. The vector r± totally 

covers the first j — I layers therefore the characteristic matrix at depth r± can be written as 

L (rx) = Lj (rx - Dj.,) ■ L(,_i) (B3) 

where 7?j_i = J^d; is the total thickness of layers up to layer j — 1 and 7v(j_i) = Lj_i • ... • Li 

1=1 

is the characteristic matrix of layers 1, .., j — 1. We note that layer j is only partially covered 
by the depth interval, which is indicated by the argument (r_|_ — Dj_i) in Eq. (B3) instead 
of the total thickness dj. It is also important to note that L (rj,) depends on the thicknesses 
and susceptibilities of all covered layers and, furthermore, that it also depends on the angle 
of grazing incidence 6. 
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Using Eqs. (2.8), (B2) and (B3) the integral (2.12) can be analytically calculated. Indeed, 
the two integrals 



result in 



with 



If = ^= / dr_L exp {-ik_i_r±) smh [k {r± - -Dj-i) FA 
V27r J 

Ij = J dr± exp {—ik'j_r±) cosh [k {r± — -Dj-i) Fj] 



aj = exp {—ik'j_Dj_i) 
Pj = exp {—ik'j_Dj_i) 



(A7)-'exp(|A7)smh(|A- 
(A7)-'exp(-|/f;)Bmh(|A: 



Here Kf — kFj ± ik'j^I. Finally the required expression has the form 



(B4a) 
(B4b) 



(B5a) 
(B5b) 



(B6a) 
(B6b) 
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■ I^-lr+r[12] .-.[22] 

^^j h ^(i-1) + h ^0-1) 



(/ + i?sp) ■ 
(B7) 

Eq. (B7) is physically the Fourier transform of the depth-profile function of the coherent 
field and we emphasize again its dependence on k^. 



APPENDIX C: DWBA OFF-SPECULAR INTENSITY FORMULAE 

As it was noted in Sections II and III, the diffuse intensity expressions (2.18) and (3.10) 
differ from the result of DWBA. Indeed, the DWBA transition matrix element {(3 \T\ a) of 
the transition between the eigenstates \a) and of the interaction- free Hamiltonian is 



{(5\T\a)^{Kl^,V24a). 



(CI) 



where T is the transition matrix, V2 is the perturbing potential, k^^ and are the re- 
tarded and the advanced solutions of the unperturbed Hamiltonian and adjoint Hamiltonian, 
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asymptotically related to \a) and respectively.^^ The interaction potential V — k'^x is 
split according to Eq. (2.4) to a sum 

where the wave equation is exactly soluble for Vi = k'^x ^-iid V2 = k'^ {x ~ x) is regarded as 
the perturbing potential. The arguments of x and x were dropped for the sake of simplicity. 
According to the concept used in this paper, the initial and final states \a) and are plane 
waves with wave vectors k and k', respectively. 

The solution of the unperturbed equation (2.3), having only the Vi potential, was given 
in Eq. (2.7). One notes that Eq. (2.7) is the retarded solution of Eq. (2.3), therefore 

< = T {k^, r^) *-exp (zk||r||) . (C3) 

The advanced solution can be given similarly, taking into account that the advanced solution 
is related to the emerging plane wave rather than to the incoming wave considered in Eq. 
(A2). Therefore using 

* (0) = + (C4) 

and Vi instead of Vi, the advanced solution can be written as 



K 
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= f (-A;l, rx) *°"*exp (ik||r||) , (C5) 



where T means that in Eq. (A5) the specular reflectivities are inverted and the adjoint 
susceptibilities are used. Applying Eqs. (CI), (C3), (C5), (3.3), (3.4), (3.8) and (3.12) the 
diffuse intensity is 

/off (E) ^k'Yl (K||) IV [r|;' {E)' rf {E) p] (C6) 
similarly to Eq. (3.10), however, the definition of being 

rf (A;^,^!,/;) = y drxf (-A;;,rx)^xr(£^)r(^±,r±)- (C7) 

different from that given in Eq. (3.11). The integral (C7) can be analytically calculated 
following the method presented in Appendix B. Comparing Eq. (C7) with Eqs. (2.12) and 
(3.11) one can reahze that the present result of Eq. (3.10) can be also obtained from the 
DWBA by taking 

f (-A;l,rx)^ Ri7exp(-iA;lrx), (C8) 
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which is vahd only for exit angles ( 9') above the critical angle. This statement is consistent 
with ignoring the scattering of the diffuse field as assumed in section II. We also note that 
approximation 

T {k^, rx) ~ / exp {ik^r^) (C9) 
together with (C8) are identical to the conventional 1®* order Born approximation. 
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FIG. 1: Simulated ^ Fe specular synchrotron Mossbauer reflectograms {6 — 29 scans) of the 
MgO/ [^^Fe (2.62nm) /Cr (1.28nm)]2Q antiferromagnetic multilayer with sublayer magnetizations 
parallel and antiparallel to the wave vector for different domain bias parameters, ij indicated in the 
figure. 
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FIG. 2: Simulated ofF-specular PNR a;— scans calculated for the 

MgO/ [5^Fe (2.62 nm) /Cr (1.28 nm)l antiferromagnetic multilayer with A = 0.4 nm and 
detector position 26 fixed at the 1/2-order antiferromagnetic Bragg peak position for various 
correlation lengths, ^ indicated in the figure. A single domain bias parameter of 77 = 0.1 was used. 
Dotted and solid lines show the BA and SDIWA simulations, respectively. 
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FIG. 3: Simulated off-specular SMR a;— scans calculated for the 

MgO/ [^^Fe (2.62 nm) /Cr (1.28 nm)]2Q antiferromagnetic multilayer with A = 0.086 nm of 
the ^^Fe Mossbauer radiation and the detector position 29 fixed at the 1/2-order antiferromagnetic 
Bragg peak position for various correlation lengths, ^ indicated in the figure. A single domain 
bias parameter of ?7 = 0.1 was used. Dotted and solid lines show the BA and SDIWA simulations, 
respectively. 
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FIG. 4: Simulated '6 

in ~ ^out' PNR diffuse intensity maps for the 
MgO/ [^"^Fe (2.62nm) /Cr (1.28nm)]2Q antiferromagnetic multilayer structure. The intensi- 
ties are normalized and shown on a logarithmic color scale. BA (a) and SDIWA (b) J"*""^ intensities 
are shown. A single domain bias parameter of 77 = 0.1 was used. 
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FIG. 5: Simulated '^in — ^out' SMR diffuse intensity maps of the 

MgO/ [^'^Fe (2.62nm) /Cr (1.28nm)]2Q antiferromagnetic multilayer structure around the 1/2- 
order antiferromagnetic Bragg peak using the BA (a) and the SDIWA (b). The intensities are 
shown on a logarithmic color scale and are normalized. A single domain bias parameter oi r] = 0.1 
was used. 
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FIG. 6: Measured (a) and simulated (b) — 2^' SMR diffuse intensity maps of the 
MgO/ p7Fe(2.62nm) /Cr(1.28 nm)l antiferromagnetic multilayer in the vicinity of the 1/2-order 
antiferromagnetic Bragg peak. The intensities are normalized and shown on a logarithmic color 
scale. Map (b) was simulated using Eqs. 3.20, 4.3 and 4.4 with ^ = 1.0 /im domain correlation 
length. The dotted white line shows the 10% level of the maximum intensity. 
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